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1.  Introduction 

In  a  recent  paper,  Lee  and  Chen  (2012,  hereafter 
LC12)  presented  numerical  simulations  of  symmetric 
and  asymmetric  hurricane  boundary  layer  structures  in 
a  fully  coupled  atmosphere-wave-ocean  model  and 
used  these  simulations  to  compare  aspects  of  the 
boundary  layer  structure  against  an  analysis  of  obser¬ 
vations.  One  of  their  main  conclusions  was  that  “the 
azimuthally  averaged  inflow  layer  tends  to  misrepresent 
the  overall  inflow  structure  in  tropical  cyclones,  espe¬ 
cially  the  asymmetric  structure”  (p.  3593).  Another 
main  conclusion  was  that  the  complicated  asymmetric 
three-dimensional  boundary  layer  structures  (attrib¬ 
uted  by  them  to  be)  due  in  part  to  the  air-sea  and  wind- 
wave  coupling  “make  it  difficult  to  parameterize  the 
atmosphere-wave-ocean  coupling  effects  without 
a  fully  coupled  model”  (p.  3593).  After  careful  exami¬ 
nation  of  their  study,  we  have  a  number  of  questions 
regarding  their  methodology,  their  interpretations  (in¬ 
cluding  their  interpretations  of  previous  literature),  and 
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their  conclusions.  Specifically,  we  inquire  about  aspects  of 
the  methodology  for  defining  the  dynamical  boundary 
layer  depth,  the  selection  of  the  boundary  layer  scheme, 
and  we  question  the  conclusions  inferred.  In  addition  to 
the  foregoing  concerns,  inaccuracies  in  their  literature 
review  are  noted  and  inconsistencies  between  their  con¬ 
clusions  and  reported  results  are  identified. 

For  many  decades,  physical  processes  across  the  air- 
sea  interface  and  within  the  atmospheric  boundary  layer 
have  been  known  to  be  essential  for  the  development 
and  maintenance  of  a  tropical  cyclone  (Ooyama  1969; 
Emanuel  1986,  1995;  Smith  et  al.  2009;  Smith  and 
Montgomery  2010;  Bryan  and  Rotunno  2009;  Bryan 
2012).  However,  the  boundary  layer  is  the  least-observed 
part  of  a  storm — in  particular,  its  turbulence  structure. 
With  the  advent  of  the  global  positioning  system  (GPS) 
dropsonde  (Hock  and  Franklin  1999),  the  mean  bound¬ 
ary  layer  structure  has  been  progressively  studied.  Pre¬ 
vious  studies  have  concentrated  mostly  on  determining 
the  boundary  layer  structure  in  an  individual  storm  (e.g., 
Kepert  2006a,b;  Montgomery  et  al.  2006;  Bell  and 
Montgomery  2008;  Barnes  2008;  Zhang  et  al.  2009;  Zhang 
2010)  with  the  hope  that  these  findings  generalize  to  other 
storms.  Recently,  Zhang  et  al.  (2011,  hereafter  Zll) 
conducted  a  composite  analysis  of  the  axisymmetric 
boundary  layer  structure  based  on  hundreds  of  GPS 
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dropsondes  from  13  Atlantic  hurricanes.  They  found  that 
there  is  a  clear  separation  between  the  boundary  layer 
depths  defined  kinematically/dynamically  and  thermody¬ 
namically.  The  kinematic  boundary  layer  depth  (i.e.,  the 
depth  of  relatively  strong  inflow  possessing  greater  than 
10%  of  the  peak  mean  inflow  velocity)  is  much  higher  than 
the  thermodynamic  height  (i.e.,  the  mixed-layer  depth  for 
virtual  potential  temperature).  They  found  also  that  the 
observed  inflow-layer  depth,  as  well  as  the  depth  of  rela¬ 
tively  strong  inflow,  tends  to  decrease  with  decreasing 
radius  toward  the  center.  These  results  support  strongly 
those  found  in  the  earlier  studies  of  specific  storms. 

The  study  by  LC12  presents  an  analysis  of  the  asym¬ 
metric  structure  of  the  atmospheric  boundary  layer  in  a 
strong  hurricane  vortex.  Based  on  comparisons  between 
simulated  and  observed  vertical  profiles  of  kinematic 
and  thermodynamic  parameters,  primarily  from  data 
collected  in  Hurricane  Frances  (2004),  they  suggested 
that  the  atmosphere-ocean  coupling  reduces  the  mixed- 
layer  depth  in  the  rear-right  quadrant  because  of  storm- 
induced  ocean  cooling  and  that  the  wind-wave  coupling 
enhances  boundary  layer  inflow  beyond  the  radius  of 
maximum  wind  speed.  They  reported  also  a  signiflcant 
front-to-back  asymmetry  in  the  depth  of  the  inflow  layer 
in  their  numerical  simulations.  LC12  emphasized  the 
complicated  nature  of  the  three-dimensional  boundary 
layer  structure  in  their  simulations.  However,  they  offered 
no  discussion  or  reasoning  why  such  asymmetric  struc¬ 
tures  would  be  important  for  hurricane  intensification  or 
mature  intensity.  We  communicate  our  questions  here. 

As  stated  above,  our  study  of  LC12  has  raised  what  we 
believe  are  substantive  questions  about  their  analysis 
method  and  scientific  interpretation  that  diminish  the 
reliability  of  some  of  their  conclusions.  Inaccuracies  in 
their  literature  review  are  noted  also. 

2.  Scientific  concerns  with  LC12’s  methodology 

Our  first  major  scientific  concern  with  LC12  is  their 
definition  of  the  hurricane  boundary  layer.  They  wrote 
“...  to  distinguish  the  inflow  layer  from  that  of  the 
mixed  layer,  here  we  refer  to  it  as  dynamic  HBL 
(DHBL)”  (p.  3577).  Thus,  they  equated  the  depth  of  the 
dynamical  hurricane  boundary  layer  (DHBL)  with  the 
height  of  zero  inflow.^  LC12  wrote  also  that  “The  bound¬ 
ary  inflow  is  a  result  of  gradient  wind  imbalance  due  to  the 
surface  friction  (Smith  1968),  and  the  top  of  the  inflow 
layer  is  defined  as  where  the  inflow  vanishes  (e.g..  Smith 
et  al.  2009;  Zhang  et  al.  2011)”  (p.  3577).  They  cited  Smith 


^  Note  that  the  definition  in  the  reply  (Lee  and  Chen  2014)  as  the 
contour  of  — 2ms~^  is  different  from  that  in  the  original  paper. 


et  al.  (2009)  and  Zll  as  if  these  two  studies  advocated 
such  a  definition  for  the  boundary  layer.  In  fact.  Smith 
et  al.  (2009)  defined  the  boundary  layer  depth  as  the 
depth  of  the  “the  shallow  layer  of  strong  inflow  near  the  sea 
surface  that  is  typically  500  m  to  1km  deep  and  arises 
largely  because  of  the  frictional  disruption  of  gradient  wind 
balance  near  the  surface”  (footnote  on  p.  1322).  This 
definition  does  not  include  the  weak  midlevel  inflow.  This 
depth  was  easy  to  distinguish  in  their  numerical  experi¬ 
ments.  In  the  same  spirit,  Zll  defined  the  boundary  layer 
height  as  the  height  at  which  the  inflow  falls  to  10%  of  the 
peak  inflow  (see  their  Fig.  5).  The  inappropriate  choice  of 
the  inflow  layer  depth,  including  the  convectively  driven 
deep  tropospheric  inflow  by  LC12,  seems  to  be  one  of  the 
reasons  for  their  misinterpretation  of  results  reported  by 
Zll  and  would  explain  why  the  DHBL  depth  shown  in 
their  Fig.  17  is  as  high  as  10  km. 

Zll  were  fully  aware  of  the  potential  limitations  of 
using  the  inflow-layer  depth  as  the  top  of  the  hurricane 
boundary  layer  and  wrote  an  entire  paragraph  articulat¬ 
ing  the  issues  involved.  On  p.  2531-2532,  they  wrote; 

Notwithstanding  the  variability  of  different  boundary 
layer  height  scales,  it  is  thought  that  the  inflow  layer  depth 
represents  the  top  of  the  hurricane  boundary  layer  better 
than  does  the  thermodynamic  boundary  layer  depth.  Di¬ 
rect  flux  measurements  in  the  outer-core  regions  of  hur¬ 
ricanes  suggest  the  turbulent  flux  transport  mainly  occurs 
in  the  inflow  layer  (Zhang  et  al.  2009).  The  budgets  and 
discussion  presented  by  Kepert  and  Wang  (2001)  and 
Kepert  (2010a)  support  the  statement  that  the  momentum 
flux  occurs  mainly  in  the  inflow  layer.  In  his  numerical 
simulations,  Kepert  (2010a)  also  showed  that  the  mo¬ 
mentum  flux  is  a  significant  part  of  the  dynamics  of  the 
layer  of  outflow  immediately  above  the  inflow  and  sug¬ 
gested  that  it  is  therefore  appropriate  to  include  at  least 
part  of  this  layer  in  the  boundary  layer. 

LC12  challenged  Zll’s  results  (see  their  last  sentence  on 
p.  3589)  without  comment  on  the  foregoing  caveats. 

A  second  major  scientific  concern  that  we  have  with 
LC12’s  study  is  a  fundamental  flaw  in  their  methodology 
of  comparing  the  model  asymmetric  structure  with  the 
observed  composite  symmetric  structure  from  many 
storms  in  Zll.  We  would  argue  that  this  comparison  is 
logically  ill  founded.  For  a  single-valued  function  of  azi¬ 
muth  0,  the  azimuthal  structure  is  a  periodic  function  of 
6  at  a  given  radius  r  and  height  z.  The  axisymmetric 
structure  is  represented  by  azimuthally  averaging  the 
field.  By  definition,  at  any  point  in  time,  the  azimuthal- 
mean  structure  is  mathematically  orthogonal  to  the 
asymmetric  structure.  For  this  reason,  Zll  never  sug¬ 
gested  (either  explicitly  or  implicitly)  that  their  azimuthal- 
mean  structure  would  represent  the  asymmetric 
structure. 
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3.  Concerns  with  LC12’s  interpretations  of  resnlts 
relating  to  past  literatnre 

Besides  their  misinterpretation  of  the  methodology  of 
defining  the  kinematic/dynamic  boundary  layer  depth 
used  by  Zll  and  Smith  et  al.  (2009),  LC12  questioned 
several  times  the  representativeness  of  Zll’s  results.  For 
instance,  they  stated  that  “Many  [storms]  have  a  much 
deeper  inflow  layer  in  parts  of  the  hurricanes  than  the 
composite  in  Zll,  which  raises  a  question  of  whether  the 
composite  inflow  can  represent  the  true  structure  in 
hurricanes”  (p.  3577).  They  go  on  to  say  that  “Unlike 
the  azimuthally  averaged  fields  shown  in  Zhang  et  al. 
(2011),  the  low-level  outflow  layer  above  DHBL  [dy¬ 
namic  hurricane  boundary  layer;  our  insertion]  exists  only 
in  the  front-left  quadrant  (Kepert  2006a)”  (p.  3578).  They 
stated  also  that  “These  features  are  different  from  the 
mean  inflow  layer  described  in  Zhang  et  al.  (2011),  but 
they  are  consistent  with  that  in  Hurricane  Georges  (1998) 
shown  in  Kepert  (2006a)”  (p.  3586).  Furthermore,  they 
stated  that  “It  also  causes  concern  regarding  the  repre¬ 
sentativeness  of  the  azimuthally  averaged  HBL  proper¬ 
ties  as  shown  in  Zhang  et  al.  (2011),  which  mask  some 
dominate  [sic]  features  in  the  inflow  depth  and  asymme¬ 
try”  (p.  3593).  These  misinterpretations  have  the  poten¬ 
tial  to  confuse  readers. 

Figure  16  of  LC12  shows  that  the  mixed-layer  depth 
decreases  with  decreasing  radius  in  all  four  quadrants. 
Moreover,  the  mixed-layer  depth  is  within  the  inflow 
layer  in  all  four  quadrants.  These  findings  are  entirely 
consistent  with  the  main  conclusions  of  Zll  and  the 
conceptual  model  proposed  by  Zll  (see  their  Fig.  12  on 
p.  2532).  Thus,  the  conclusion  made  by  LC12  in  terms  of 
the  behavior  of  the  mixed-layer  depth  as  compared  to 
Zll  is  not  supported  by  their  own  results.  Furthermore, 
their  Fig.  16  shows  that  the  simulated  strong  inflow  layer 
is  generally  well  below  1-2  km  for  all  four  quadrants 
relative  to  the  storm  motion,  again  supporting  the  con¬ 
ceptual  model  of  Zll.  Note  that  Zhang  et  al.  (2013) 
conducted  analyses  of  the  asymmetric  boundary  layer 
structure  relative  to  the  environmental  shear,  sup¬ 
porting  also  Zll’s  conceptual  model  for  boundary 
layer  height  variations,  even  considering  variations  in 
azimuth. 


4.  Concerns  with  LClZ’s  conclnsions 

Finally,  we  question  the  conclusion  drawn  by  LC12  in 
terms  of  the  fully  coupled  model  in  the  last  sentence  of 
the  abstract  (and  echoed  in  their  conclusions)  that  “the 
complex,  three-dimensional  asymmetric  structure  in  both 
thermodynamic  and  dynamic  properties  of  the  HBL  in¬ 
dicates  that  it  would  be  difficult  to  parameterize  the 


effects  of  air-sea  coupling  without  a  fully  coupled  model” 
(p.  3576).  Previous  studies  have  shown  very  encouraging 
comparisons  between  model  and  observed  boundary 
layer  mean  and  turbulent  structure  without  using  a  fully 
coupled  model  (e.g.,  Nolan  et  al.  2009a,b;  Gopalakrishnan 
et  al.  2013).  Also,  in  an  idealized  study  examining  both 
a  stationary  and  a  moving  vortex,  Thomsen  et  al.  (2014) 
showed  that  tropical  cyclone  intensification  and  mature 
intensity  are  not  sensitive  to  stochastic  variations  in  the 
drag  coefficient  whose  variations  lie  within  a  reason¬ 
ably  wide  envelope.  Finally,  in  a  recent  study,  Cione 
et  al.  (2013)  found  evidence  to  suggest  that  processes  in 
the  atmospheric  boundary  layer  are  more  responsible 
for  intensity  change  than  the  ocean  response  based  on 
an  analysis  of  a  large  number  of  buoy  observations. 
Together,  the  foregoing  results  suggest  that  wind-wave 
coupling  may  not  be  necessary  for  simulating  hurricane 
intensity. 

5.  Discussion  and  recommendations 
regarding  LC12 

It  should  be  noted  that  most  of  LC12’s  conclusions 
were  drawn  using  one  planetary  boundary  layer  (PBL) 
scheme  (the  Blackadar  scheme)  in  a  comparison  of  three 
deterministic  calculations  with  distinct  model  configura¬ 
tions  of  Hurricane  Frances  (2004).  These  configurations 
included  an  atmospheric  model  only,  an  atmosphere  and 
ocean  model,  and  an  atmosphere-wave-ocean  model.  As 
pointed  out  by  Braun  and  Tao  (2000)  and  more  exten¬ 
sively  in  a  recent  study  by  Smith  and  Thomsen  (2010), 
detailed  aspects  of  the  simulated  hurricane  boundary 
layer  structure  are  generally  sensitive  to  the  selection  of 
PBL  scheme.  We  note  also  that  the  Blackadar  scheme 
used  a  constant  vertical  mixing  length,  which  is  not  re¬ 
alistic  according  to  direct  flux  measurements  in  the  hur¬ 
ricane  boundary  layer  by  Zhang  and  Drennan  (2012). 

We  have  an  additional  concern  with  LC12’s  method¬ 
ology  of  using  a  single  deterministic  calculation  for  each 
model  configuration  in  light  of  the  variability  associated 
with  the  stochastic  nature  of  deep  convection.  For  ex¬ 
ample,  it  has  been  shown  in  a  recent  paper  by  Thomsen 
et  al.  (2013,  manuscript  submitted  to  Mon.  Wea.  Rev.) 
that  this  variability  may  lead  to  substantial  fluctuations 
in  the  asymmetric  inflow  structure  on  time  scales  as 
short  as  15  min.  Thus,  without  suitable  time  averaging  of 
the  model  or  observational  data,  it  is  unclear  whether 
the  findings  of  LC12  regarding  the  asymmetries  are 
robust. 

Given  all  of  the  issues  discussed  here,  further  analysis 
and  scientific  interpretation  of  existing  numerical  hur¬ 
ricane  simulations  following  the  same  methodology  for 
defining  the  boundary  layer  height  scales  as  used  in  the 


July  2014 


CORRESPONDENCE 


2785 


recently  published  boundary  layer  studies  cited  above  is 
strongly  encouraged. 
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